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Abstract: Functional structure and diversity of soil free-living nematodes in a desert environment depend 
on plant gender and sampling site. The objective of this study was to compare the composition, 
abundance and tropic group of soil free-living nematodes in the upper 0—10 cm soil layer under the male 
and female Acanthosicyos horridus Welw. ex Hook. f. plants and in the inter-shrub open areas (control) in the 
Namib Desert, Namibia in April 2015. Soil moisture, organic matter (OM) and pH was also analyzed. 
Free-living nematodes were extracted from 100 g soil using the Baermann funnel procedure, and total 
number of nematodes was counted under a microscope. Community composition and diversity of soil 
free-living nematodes were analyzed using 188 rDNA sequences. Results indicated that a total of 67 
groups, including 64 species, 2 genera and 1 family were identified. Feeding behavior of 58 species were 
identified as follows: 15 bacteria-feeding species, 12 fungi-feeding species, 10 plant-parasite species, 5 
omnivorous-predator species, 8 animal-parasite species, 5 invertebrate-parasite species and 3 non-free-living 
nematodes, known as marine species. Moreover, soil free-living nematodes were found to be affected by 
sampling locations and plant gender, and community composition and density of these nematodes were 
strongly influenced by soil OM content. Result confirmed that spatial location and plant cover were main 
factors influencing the diversity of soil free-living nematodes. Moreover, molecular tools were found to be 
very useful in defining the richness of soil non-free-living nematodes. In conclusion, the results elucidated 
the importance of biotic variables in determining the composition and abundance of soil free-living 
nematodes in the Namib Desert, Namibia. 
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1 Introduction 


The common feature unique to all deserts is water availability that influences substrate 
composition and structure (Whitford, 2002). In the desert ecosystem with a high temperature, soil 
moisture (SM) is unpredictable in time, space and amount. SM is critical for the biotic community, 
and enables it to fulfill its biological functions via adaptation to such arid environment (Crawford, 
1981; Evenari et al., 1982; Wardle, 2002; Whitford, 2002). 

In the Namib Desert, one of the largest sandy deserts in the world (Laity, 2008), edaphic factors 


result in the aridity. Despite the relatively low amount of precipitation («200 mm), soil substrate 
(deep layer of sandy soil) enhances water filtration and results in the lack of water availability in 
the upper soil layer. As a result, moisture level in the upper soil layer combined with the extreme 
fluctuation in temperature drive soil biota to develop an eco-physiological adaptation due to the 
unpredictability of food availability (Evenari et al., 1982; Steinberger et al., 1989; Whitford, 
2002). Moreover, the Namib desert is also known as one of the six coastal fog deserts whose 
water balance is based on the inputs of fog and dew, in addition to the precipitation (Eckardt et al., 
2013). These two additional sources of water play important roles in the above-ground biotic 
community (Seely et al., 2005) and are of great importance to the below-ground biotic community 
(soil microflora, microfauna and mesofauna) (Rodriguez-Zaragoza et al., 2005; Stomeo et al., 
2012; Frossard et al., 2015). However, little is known about soil free-living nematodes and the 
mesofauna-plant cover interactions in this desert system. 

Soil free-living nematodes are known to constitute the most abundant and diverse group of 
multicellular animals among soil biota (Lambshead, 2004). The phylum of soil free-living 
nematodes is estimated to contain over one million species inhabiting aquatic and terrestrial 
ecosystems. These nematodes play a wide range of ecosystem function and service through 
interactions with the other groups of organisms in soils or sediments. Moreover, their interactions 
with other biota are not limited to feeding on bacteria, fungi, plant cells and omnivores. They also 
feed on other nematodes, heterotrophic protists and arthropod hosts (Yeates et al., 1993; Yeates 
and King, 1997; Levi et al., 2012). The activity, density and diversity of soil free-living 
nematodes in natural environments were found to be strongly influenced by soil 
physical-chemical characteristics and food resources (Pen-Mouratov et al., 2004; Fitoussi et al., 
2016). The abundance and diversity of bacteria-feeding and fungi-feeding nematodes in a desert 
are the basis of distribution of plant litter (Freckman and Mankau, 1977; Steinberger et al., 1988). 
Therefore, it is necessary to study the density and diversity of soil free-living nematodes of 
common plants in a desert system. 

Acanthosicyos horridus Welw. ex Hook. f. is widely distributed in the Namib Desert (Berry, 
2003). The male and female plants of A. horridus are equally distributed (Cloudsley-Thompson, 
1996). Moreover, the female plants of A. horridus produces succulent melon fruits that might 
contribute the accumulation of organic matter in its microhabitat. The objective of this study was 
to examine the effects of male and female A. horridus plants in the Namib Desert on the 
composition, abundance and tropic group of soil free-living nematodes. We hypothesized that 
abiotic components, e.g., SM, organic matter (OM), soil pH and plant gender play important roles 
in the composition, abundance and trophic group of soil free-living nematodes. 


2 Study area and methods 
2. Study area 


The Namib Desert is hyper-arid, which is known as a detritus-driven ecosystem in Namibia, 
Africa. The daily mean temperature ranges from 15?C to 55?C (Seely and Louw, 1980). The mean 
annual precipitation is 14.1 mm. In addition, the annual input of fog equals to 69.2 mm. The 
average number of precipitation is 59 d, ranging from 45 to 138 d (Gobabeb Research and 
Training Center (GRTC) First Order Meteorological Station, Namibia). We collected soil samples 
in the early morning in order to avoid the heating and drying of dew. The vegetation in the Namib 
Desert is covered with sparsely grasses, such as Stipagrosti sabulicola and Stipagrostis 
gonastachys, Arthraerua leubnitziae, Zygophyllum stapfii and A. horridus (Louw and Seely, 
1982). The perennial A. horridus plants that are totally leafless, attain a height of about 1.5 m, and 
have deep roots that are able to reach underground water. The female plants produce spiny fruits 
in February and April, which are eaten by animals and collected by local residents as an important 
food source. 


2.2 Soil sample collection 


Soil samples were collected from the following three sites along a 150-km distance: (1) Kuiseb 


Delta (hereinafter referred to as Delta) site (23°06.968'N, 14?28.640'E; 26 m a.s.l.) is located in a 
green belt (1.5 km from the seashore of the Atlantic Ocean). The mean amount of precipitation is 
99 mm; (2) Gobabeb site (23?05.903'N, 15?03.109'E; 418 m a.s.l.) is located in the vicinity of the 
Gobabeb Desert Research Station, about 60.0 km from the seashore of the Atlantic Ocean. The 
mean amount of precipitation ranges from 100 to 199 mm; and (3) Far East (FE) sand dunes 
(23?40.315'N, 15°39.251'E; 872 m a.s.l.) is located at a distance of 150.0 km from the seashore of 
the Atlantic Ocean, with the mean amount of precipitation ranging from 200 to 300 mm. 


2.3 Soilanalysis 


Three sampling sites are located at distances of 1.5, 60.0, and 150.0 km from the seashore of the 
Atlantic Ocean. Soil samples at each site were collected from the upper 0—10 cm (n=4) layer 
beneath the male and female A. horridus plants, and in the inter-shrub open areas as a control 
(total of 36 samples) in April 2015. Each soil sample was placed in an individual plastic bag and 
took back to laboratory. 

SM (46) content was gravimetrically determined by drying soil samples for 72 h at 105°C. OM 
(96) content was determined by igniting samples at 490°C for 8 h. Soil pH was measured with a 
glass electrode using a 1:2 soil:water ratio. Bulk density (BD, g/cm?) was determined as the 
weight of soil in a given volume. For the determination of soil water-holding capacity (WHC, 46), 
100 g soil samples were flooded with tap water in a bottom-perforated vessel for 5 min. 

Soil free-living nematodes were extracted using the Baermann funnel procedure (Cairns, 1960), 
and total number of nematodes was counted under a microscope. The nematodes were transferred 
to an Eppendorf tube after being centrifuged for 10 min at 10,000 rpm to reduce the amount of 
remaining water. The obtained pellets were placed at —20?C until identification using molecular 
tools. DNA was extracted from the pellets using the PowerSoil DNA Isolation Kit (MO-BIO 
Laboratories, Inc., Carlsbad, USA). The eluted DNA was stored at —20?C until it was used as a 
polymerase-chain-reaction (PCR) template. The 3' end of the 18S rDNA was amplified using 
primer pair NF1 (5'-GGT GGT GCA TGG CCG TTC TTA GTT-3) and 18Sr2b (5-TAC AAA 
GGG CAG GGA CGT AAT-3) (Porazinska et al., 2009). Reactions were carried out at a 25.00-uL 
final volume including 15.75 uL ultrapure water, 5.00 uL buffer, 0.50 uL dNTPs, 1.25 uL F 
primer, 1.25 uL R primer, 0.25 uL Phusion High Fidelity DNA Polymerase (100 U, 2 U/uL, 
Thermo Scientific, USA) and 1.00 uL DNA from each sample. Cycling was performed by Veriti 
96-Well Thermal Cycler (Applied Biosystems, Carlsbad, California, USA). The PCR program 
involved a denaturation step of 30 s at 98?C, a denaturation step of 10 s at 98?C, an annealing step 
of 30 s at 56?C, an extension step of 1.5 min at 72?C and a final extension step of 10.0 min at 72?C. 
Steps 2—4 were repeated for 20 cycles. Each PCR product was diluted with ultrapure water 1:100, 
and underwent both vortex spinning and a second PCR amplification using F primer, with 
barcodes for each sample. A negative control was included in each amplification. 

For 18S rDNA sequences, we carried out amplification using the forward primer barcode 
sequence NTF1 IT1-NTF1 IT37: 5'-CCA TCT CAT CCC TGC GTG TCT CCG ACT CAG CTA 
AGG TAA CGG TGG TGC ATG GCC GT-3', and the reverse primer sequence 18Sr2bIT: 5'-CCT 
CTC TAT GGG CAG TCG GTG ATT ACA AAG GGC AGG GAC GTA AT-3'. The PCR program 
involved a denaturation step of 3 min at 98°C, a denaturation step of 30 s at 98°C, an annealing 
step of 30 s at 55°C, an extension step of 1 min at 72°C and a final extension step of 5 min at 
72°C. A 5 uL portion of all PCR products, along with 1 kb DNA ladder (Thermo Scientific, USA) 
was electrophoresed in 1% agarose gel (SeaKem, Lonza, Rockland, ME, USA) and stained with 1 
uL GelStar. No PCR products were observed in the negative control, and the correct size of the 
positive band was observed. The remaining 20 uL DNA of each sample were electrophoresed in 
0.5% low-melt agarose (CSL-AG100), and stained with a 3-uL gel. The GeneJET Kit (Thermo 
Scientific, USA) was used for the recovery and purification of DNA from low-melting agarose. 
After cutting the desired band length of DNA (400 bp), all DNA sample products underwent ion 
torrent sequencing of 18S rDNA (Fitoussi et al., 2016). 

We processed the recovered sequences of DNA in QIIME (Quantitative Insights into Microbial 


Ecology) version 1.7.0 (Caporaso et al., 2010a, b) using the pipeline for analyzing 188 rDNA 
sequence data. Default settings, along with UClust algorithm, were used for de novo selection of 
operational taxonomic units (OTUs) at 97% similarity. The SILVA 111 release was used as a 
reference for taxonomic assignments of OTUs (Yilmaz et al., 2014). Since the SILVA 111 
database contained only an essential number of nematode reads, we sought to reclassify the reads 
assigned to Nematoda in the QIIME classification by the Basic Local Alignment Search Tool 
(BLAST) similarity analysis using a custom-made nematode 18S rDNA sequences. The database 
enabled us to obtain a better insight into the community of nematodes. The R/Bioconductor 
Package Phyloseq version 1.16.2 was used to draw phylogenetic tree from the QIIME data. All 
measurements were assessed using data for OTUs with at least 5 reads in the overall dataset. 


2.4 Data analysis 


All data were subjected to statistical analysis of variance (ANOVA) using the SAS model 
(Duncan's multiple test and Pearson's coefficient), and were used to evaluate differences between 
separate means. ANOVA was followed by Tukey's HSD test to establish the significance of 
differences between plot areas using the statistical package Statistica 4.3 (StatSoft Inc., Tulsa, OK, 
USA). Differences obtained at level of P<0.05 were considered significant. 

We used RDA to correlate the composition of soil free-living nematodes with environmental 
variables. Partial RDA were used to determine the effect of variable as covariable, that of soil 
water content with the remaining variables as covariables, that of WHC with the remaining 
variables as covariables, and that of conductivity (cond) with the remaining variables as 
covariables, etc. (Liu et al., 2012). 

The data and Monte Carlo reduced model tests with 499 unrestricted permutations were used to 
statistically evaluate significance of the first canonical axis and all canonical axes combined. To 
meet requirements of the Monte Carlo reduced model test, we square-root-transformed the data of 
taxonomic group. DCA, RDA and partial RDA were carried out using CANOCO software for 
Windows 4.5 (Microcomputer Power, Ithaca, NY, USA). 


3 Results 


3.1 Abiotic variables 


Mean values of different soil physical parameters at three sites and three locations, i.e., male and 
female A. horridus plants and the inter-shrub open areas as the control at each site are presented 
in Table 1. 

SM contents were low and ranged from 0.24% at Delta to 0.49% at FE sand dunes. Significant 
difference in SM (P«0.05) content at FE sand dunes was found. SM content in three locations at 
each site did not significantly differ (Table 1). 

OM contents ranged from 0.26% at FE sand dunes to 0.52% at Gobabeb site. OM content was 
significant (P«0.05) higher at FE sand dunes than at the other two sites. No significant differences 
were found at the locations of Gobabeb and FE sand dunes, while at Delta, OM content beneath 
the male and female plants was significant (P«0.05) higher than that of control. BD was 
significant (P«0.05) higher at Gobabeb compared with the other two sites, while WHC was 
significant (P«0.05) lower at Gobabeb. An attempt was made to determine the effects of 
male/female plants and control (i.e., three locations) at each site on soil physical-chemical 
characteristics (Table 1). For BD, significant (P«0.05) differences, such as control»female» male, 
were found only at Delta. As for WHC, no significant differences were found at three locations of 
FE sand dunes, while WHC beneath the male A. horridus plants was significant (P«0.05) higher 
than those of beneath the female A. horridus plants and control. 


Table 1 ANOVA results for soil physical-chemical characteristics at different sites and locations 
Site Location SM (96) OM (96) BD (g/cm?) WHC (96) 


Delta Cont. 0.29+0.06° 0.220.05” 1.54+0.02* 0.22+0.005° 


Fem. 0.21+0.10* 0.35+0.07* 1.38+0.02° 0.26+0.01° 


Male 0.22:40.22° 0.44+0.03" 1.31+0.03° 0.28:+0.018 
Mean 0.24+0.13° 0.330.108 1.4140.108 0.25+0.02" 

Cont. 0.2430.04* 0.53+0.02" 1.51+0.02" 0.21+0.005° 

Gobabeb Fem. 0.21+40.09" 0.45+0.04" 1.5140.04" 0.230.019 
Male 0.3840.14* 0.65+0.05" 1.49+0.02" 0.240.018 

Mean 0.26+0.108 0.52::0.08^ 1.46+0.144 0.220.018 

Far East (FE) Cont. 0.530.128 0.26+0.18" 1.4140.03" 0.26+0.018 
sand dunes Fem. 0.4740.11° 0.2340.04" 1.4040.02* 0.2640.005" 
Male 0.46+0.05° 0.29+0.02" 1.36+0.03" 0.2740.01* 

Mean 0.49+0.09^ 0.260.108 1.39+0.038 0.26+0.014 


Note: Fem., female A. horridus plants; Cont., control (inter-shrub open areas); SM, soil moisture; OM, organic matter; BD, bulk density; 
WHC, water-holding capacity. Different lowercase letters within a column indicate significance among different locations at P<0.05 
level. Different uppercase letters within a column indicate significance among different sites at P<0.05 level. Mean+SD. 


3.2 Total numbers and tropic group of soil free-living nematodes 


Density of soil free-living nematodes at different sites fluctuated from 0 to 1585 individuals/100 g 
dry soil (Fig. 1). The maximal density of soil free-living nematodes was found at FE sand dunes, 
i.e., 1585 individuals/100 g dry soil in the sample collected beneath the male plants. At Delta, 
densities of soil free-living nematodes were low, reaching 119 and 124 individuals/100 g dry soil 
in the samples beneath the female and male A. horridus plants, respectively (Fig. 1), which was 
significant lower than those of the other two sites (P<0.001). Densities of soil free-living 
nematodes were not significant at control among three sites. Similarly, no significant difference 
was found between male and female A. horridus plants, demonstrating no significant relationship 
between soil free-living nematodes abundance and plant gender. 


Density of soil free-living nematodes 
(individuals/100 g dry soil) 


Cont. Fem. Male} Cont. Fem. Male} Cont. Fem. Male 
Delta Gobbabeb FE sand dunes 


Fig. 1 Density of soil free-living nematodes at different sites (Delta, Gobabeb and Far East (FE) sand dunes) 
and locations. Fem., female A. horridus plants; Cont., control (inter-shrub open areas). Different lowercase letters 
indicate significant differences among three locations at P«0.05 level. 


A total of 6915 quality-filtered sequences were obtained and clustered into 348 OTUs at 97% 
similarity. We calculated relative abundances of metazoa, fungi, algae and protista groups based 
on the 348 OTUs. At control location of FE sand dunes, the highest relative abundance was 
Protista (Fig. 2). The Protista kingdom was represented by Alveolata, Discoba and Rhizaria phyla, 
in which Discoba is known to be adapted to saline habitats. Algae, known to be a diverse group of 
photosynthetic organisms, represented by Chloroplastida and Stramenopiles-Chrysophyta at a 
relatively low abundance of 0.3% to a maximum abundance of 23.0% at Delta beneath the male A. 
horridus plants. Fungi were represented by Ascomycota and Basidiomycota, with relatively low 
abundances ranging from 0.2% to 46.0% at Gobabeb. In general, a significant variation among 


different eukaryotic organisms was found at different sites and locations (Fig. 2). 

We found a total of 236 soil free-living nematodes based on the 348 OTUs. Low abundance 
OTUs (those containing less than 0.05% of total reads) were filtered out. Due to broad 
distribution of sequences among taxa, we selected representative sequence at a 95% similarity 
level from each OUT (Table 2). The overlap percentage of soil free-living nematodes among 
gender and control was shown in Table 2. As a result, a total of 67 groups, including 64 species, 2 
genera and 1 family were identified (Table 3). 
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Fig. 2 Relative abundances of metazoa, fungi, algae and protista groups at different sites and locations based on 
18S rDNA sequences. Fem., female A. horridus plants; Cont., control (inter-shrub open areas). 


Table 2 Percentages of soil free-living nematodes at different sites and locations, and overlap percentage among 
plant gender and control locations 


Location Delta (96) Gobabeb (96) Far East sand dunes (96) 
Male 17 50 0 
Female 33 0 13 
Control 33 0 37 
Male and female 0 0 37 
Male and control 0 25 0 
Female and control 0 13 13 
Male, female and control 17 13 0 


Unique species of soil free-living nematodes at Delta site was Ditylenchus destructor that was 
only found beneath the male A. horridus plants. The total number of soil free-living nematodes 
that was only occurred at Gobabeb was 46, of which 2 were found at control, 5 beneath the 
female A. horridus plants and 39 beneath the male A. horridus plants. Only 1 species, i.e., 
Ditylenchus sp. 5JH-2014 that was only occurred beneath the female A. horridus plants at FE 
sand dunes. The highest number of soil free-living nematodes was found beneath the male A. 
horridus plants at Gobabeb (64), followed by beneath the male A. horridus plants at FE sand 
dunes (24) and male A. horridus plants at Delta (17). No significant differences were found 
beneath the female A. horridus plants among three sites. A relatively low number of soil 
free-living nematodes were found at control, i.e., 14 at Delta, 22 at Gobabeb and 10 at FE sand 
dunes. 

Among the 64 species of soil free-living nematodes, we found 18 representatives (Table 4) of 
different habitats, such as mammalian parasites, avian parasites, microphages, marine-living 
systems, pathogens and parasitic nematodes. 

Table 3  Taxa and trophic groups of soil free-living nematodes using 18S rDNA sequences and BLAST 
similarity analysis at different sites and locations 

Delta Gobabeb Far East sand dunes — Trophic 

Cont. Fem. Male Cont. Fem. Male Cont. Fem. Male group 


Nematode species 


Ditylenchus destructor - E t FF 
Number of species - - 1 
Acrobeles ciliatus + - + BF 
Amblydorylaimus isokaryon = = + OP 
Anisakis simplex = = +++ AP 
Aphelenchoides (genus) = = + FF 
Aphelenchoides sp. USO2 = + + FF 
Brugia timori = = +++ ™ 
Bursaphelenchus abruptus - s t FF 
Caenorhabditis elegans E - + BF 
Cephalobus cubaensis = = + BF 
Cooperia oncophora E = + AP 
Echinocephalus overstreeti = = ++ 
Enchodelus longispiculus i = + PF 
Enterobius vermicularis = = + AP 
Fergusobia sp. 281 = = + FF 
Fergusobia sp. 469 i x + FF 
Halalaimus sp. TCR26 = = + MS 
Halalaimus sp. TCR93 - - + MS 
Isolaimium multistriatum - s + MS 
Krefftascaris sharpiloi - z t AP 
Longidorus jonesi = + +++ PP 
Meloidogyne javanica a = + PP 
Meloidogyne sp. Pak.P.R.m2 = = PP 
Meloidogyne sp. VO-2014 = = PP 
Panagrolaimidae sp. NK-2011b - - BF 
Parastrongyloides trichosuri = + +++ OP, AP 
Passalurus ambiguus - = +++ OP 
Pseudhalenchus minutus = = + 
Rhyssocolpus paradoxus = + OP 
Rhyssocolpus vinciguerrae = = + OP 
Robustodorus megadorus = E t 
Romanomermis culicivorax + + + 
Schistonchus caprifici = = + PP 
Steinernema feltiae E = + PI 
Steinernema sp. T51 = = + PI 
Syphacia muris E = + 
Tricoma sp. 2 CYC-2009 = = + 
Trophotylenchulus sp. TSH-2009 = = + 
Xiphinema krugi = = ++ 
Xiphinema setariae - = +++ PP 
Number of species 2 5 39 PP 
Ditylenchus sp. 5 JH-2014 = + = FF 
Number of species 0 1 0 
Acrobeles sp. = + + + = + = + + BF 
To be continued 
Continued 
] Delta Gobabeb Far East sand dunes — Trophic 
Nematode species 
Cont Fem. Male Cont. Fem. Male Cont. Fem. Male group 
Acrobeles sp. MA-2012 + + + + + + + = + BF 


Acrobeloides cf. buetschlii 1 JH-2012 +++ +++ - + j ++ t + + BF 
Acrobeloides maximus +++ +++ + +++ s +++ - t ++ BF 
Aphelenchus avenae - + 5 7 = + - + ++ FF 
Aporcelaimellus sp. JH-2004 + + + ++ +++ + + + 
Bursaphelenchus anatolius = + & = + ++ - + FF 
Bursaphelenchus penai + = + + + = + FF 
Cephalobidae (family) + ++ E +++ +++ +++ + +++ ++ BF 
Cylicostephanus goldi + = = + + + +++ - 4 AP 
Deladenus proximus - + $ + + +++ - + + PW 
Haemonchus placei ++ + +++ +++ +++ +++ + pa PPR 
Hexamermis albicans + + + + + 5 = EA 
Hexatylus sp. 1 JH-2014 = + - + + - + + FF 
Malenchus pressulus = + - + + - + + PPR 
Meloidogyne incognita - ? + 7 - + - + fe PP 
Panagrolaimus paetzoldi = i + = - ++ - + T BF 
Panagrolaimus sp. ASO1 * +++ +++ ++ +++ +++ + +++ +++ BF 
Panagrolaimus sp. AS03 + +++ +++ - = + - + + BF 
Panagrolaimus sp. SN103 z + ++ + + + - + + BF 
Paratylenchus goldeni PP +++ +++ + +++ + +++ + ++ ++ PP 
Paratylenchus labiosus - + Š - + - + ES PP 
Turbatrix aceti 7 - + + + + à + + AP 
Tylenchulus semipenetrans + ++ + + + + = + + FF 
Uncultured Rhabdolaimus (genus) + z + ++ + +++ +++ + + BF 
Number of species 14 20 16 20 18 25 10 22 24 
Total number of species 14 20 17 22 23 64 10 23 24 


Note: BF, bacteria-feeding; FF, fungi-feeding; PP, plant-parasite; AP, animal-parasite; PW, parasites on woodwasp; PPR, 
parasite/pathogenic of ruminants; EA, endoparasites in arthropods; MS, marine species; TM, transmitted by mosquitoes; PI, parasites on 
invertebrates; PF, feeds on mosses and lichens; Cont., control; Fem., female. +, appears in one replica; ++, appears in two replicas; +++, 
appears in all samples; -, no occurrence. 


Of the 64 species of soil free-living nematodes, we identified feeding behavior of 58 species, 
i.e., 15 bacteria-feeding (BF) species, 12 fungi-feeding (FF) species, 10 plant-parasite (FP) 
species, 5 omnivorous-predator species, 8 animal-parasite (AP) species, 5 invertebrate-parasite 
(PI) species and 3 non-free-living nematodes, known as marine species (MS). One species, 
Enchodelus longi spicules that feeds on mosses and lichens (PF), might belong to plant feeding 
trophic group (Table 3). 


3.3 Relative abundance of genera of soil free-living nematodes and Venn diagram 


The results of relative abundance of genera of soil free-living nematodes (Fig. 3) showed a broad 
diversity, especially the relative abundance of Panroglimus genera that showed a higher 
percentage beneath the male and female plants at three sites. 

A Venn diagram can be used to illustrate similarity and dissimilarity of soil free-living 
nematodes at each locations and possible relationships in plant gender. In three sets of circles for 
each gender and control at each site, we can see intersection or overlap of the three circles (Fig. 4). 
At Gobabeb and FE sand dunes, we found a total of 8 groups of soil free-living nematodes 
(10096), while only 6 groups were found at Delta (Fig. 4). This result could be attributed to the 
effect of abiotic stress, i.e., high salinity in the soil of Delta (data not shown). 

Table 4 Family and habitat of these species of soil free-living nematodes in this study 


Species/genus Family Habitat Reference 


Alaninema Aphrophridae Intestinal parasites in slugs Ivanova et al. (2013) 


Amblydorylamus Dorylaimidae Nematode from Antarctica Elshishka et al. (2015) 
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Brugia-B. malayi Filariidae Present only in Southeast Asia Triteeraprapab et al. (2001) 
Bursaphelenchus Parasitaphelenchidae Obligate mycrophages Ryss et al. (2005) 
Cooperia Cooperiidae Intestinal parasites in cattle Dorny et al. (1997) 
Cylicostephanus—C. goldi Strongylidae Intestinal parasites of horses Bucknell et al. (1995) 
Deladnus Neotylenchidae Parasitic nematodes Bedding (2009) 
Deladenus siricidicola Neotylenchidae Biological control agent of woodwasp Bedding (2009) 
Dicelis Drilonematoidea Parasitic earthworm Spiridonov (1992) 
Drasico nemoralis Endemic earthworm from far East Russia Ivanova et al. (2014) 
Enerobius Oxyuridae Intestinal parasite in humans Brown (2006) 
Epsilonema Epsilonematidae Marine species der n CE 
Eubosrichus Desmodoridae Marine species Polz et al. (1999) 
Gongylongema Gongylonematidae Bird and mammal parasites Soulsby (1982) 
Greeffiella Desmoscolecoidea Free-living marine nematode AMOR poo 
Haemonchus Trichostronylidae Pathogenic nematodes of ruminants Fleming et al. (2006) 
Halalaimus Oxystominidae Marine species Turpeenniemi (1998) 
Halomonhystera Monhysteridae Marine species Tchesunov et al. (2015) 
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Fig. 3 Relative abundance of genera of soil free-living nematodes using 18S rDNA sequences and BLAST 
similarity analysis at three sites and locations 


The Venn diagram also showed that at Delta, 5096 species of soil free-living nematodes 
occurred beneath the male and female A. horridus plants, while 3396 species occurred in control, 
the remaining 17% species occurred in the union of M N FNC (M, male; F, female; C, control; Fig. 
4). At Gobabeb, 50% species occurred beneath the male A. horridus plants and the remaining 
50% species occurred in the unions of MNC, FNC and MO EC. At FE, only 13% species were 
found beneath the female A. horridus plants, and 37% species at control. The remaining 50% 
species were found in the unions of MMF and FNC. 
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Fig. 4 Venn diagram showing similarity and dissimilarity of soil free-living nematodes at each locations and 
possible relationship in plant gender. Number of 0—4 means groups of soil free-living nematodes at the phylum 
level. 


3.4 Relative contributions of environmental factors to taxonomical composition 


Result of RDA showed that five environmental variables accounted for 33.0% of total variation 
(Table 5). Separate partial RDA showed that environmental variables differed in their influences 
on the composition of nematode community (Table 5; Figs. 5 and 6). 

In the five environmental variables, only the effect of OM was significant (P«0.05). Axis 1 
(F=7.05, P=0.012) showed a clear environmental gradient from Delta to FE sand dunes, which 
separated the community between Delta and FE sand dunes (Fig. 5). Axis 2 showed a clear 
environmental gradient from Gobabeb to Delta and FE sand dunes, which separated the 
community between the sites of Gobabeb and Delta+FE sand dunes. 

In the five environmental variables, only the effect of OM was significant (P«0.05). Axis 1 
(F=7.05, P=0.012) showed a clear environmental gradient from control to male, thus separating 
the community between control and male locations (Fig. 6). 


Table 5 Redundancy analysis (RDA) of environmental variables based on sampling site and plant gender 


Type Variable Initial conditional effect MCR (96) F P 
Sampling site OM 0.09 9 2.78 0.034* 
WHC 0.08 8 2.23 0.086 
Cond 0.08 8 2.23 0.078 
pH 0.07 7 2.20 0.114 
SM 0.01 1 0.35 0.824 
Total 33 
Plant gender OM 0.09 9 2.78 0.028" 
WHC 0.08 8 2.23 0.092 
Cond 0.08 8 2.23 0.092 
pH 0.07 7 2.20 0.104 
SM 0.01 1 0.35 0.770 
Total 33 


Note: OM, organic matter; WHC, water-holding capacity; Cond, conductivity; SM, soil moisture; MCR, multivariate correlation ratio. ig 


indicates significance at P«0.05 level. 


4 Discussion 


Although soil free-living nematodes are known to be one of the most diverse metazoans in 
extreme environments, they have been poorly investigated, especially in hot deserts. The ability to 
use molecular tools allow us to spend less time in identifying soil free-living nematodes present in 
soils, and to improve our understanding of soil nematodes composition, dynamic and dispersion 
in extreme dry environments (Freckman et al., 1987; Barrett et al., 2005; Treonis and Wall, 2005). 
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Fig. 5 Redundancy analysis (RDA) elucidating effects of sampling sites (D, Delta; E, Far East sand dunes; G, 
Gobabeb) on soil free-living nematodes. The number 1—23 means the replicates for sampling sites. 
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Fig. 6 Redundancy analysis (RDA) elucidating effects of plant gender and control (F, female; M, male; C, 
control) on soil free-living nematodes. The number 1—23 means the replicates for plant gender and control. 


In the present study, we assessed the effect plant gender of A. horridus on soil free-living 
nematodes along a 150.0 km distance from the Atlantic Ocean to inland in the Namib Desert 
(Yeates, 2003). Sparse vegetation covered in the Namib Desert created a "fertile island" that acted 
as food-water shoppers for both the below-ground community and the above-ground biota. In 
general, abiotic components including SM, OM, BD and soil WHC are low in xeric desert 
ecosystems mainly due to the sandy soil substrate (Yu et al., 2016). Therefore, primary production 
of vegetation and soil properties will determine the consumers, for example, soil biota that 
obtained food resources from them, and soil biotic diversity is dependent on these food resources 
(Yeates, 1971; Coleman and Crossley, 1996; Whitford, 2002; Adl, 2003). 

A significant difference in the total number of soil free-living nematodes was found along the 
150.0 km distance, with a strong influence triggered by plants and SM, as presented by 
Pen-Mouratov et al. (2010) and Fitoussi et al. (2016). Bacterial-feeding soil free-living nematodes 
were found to comprise over 50% of the samples. Their short generation time along with high 
metabolic activities and relatively high reproduction rates, allowed them to fulfill their biological 
function in an extreme xeric environment by exploiting a short period of SM availability 


(Pen-Mouratov et al., 2010). These findings were similar to those of Freckman et al. (1987) and 
Steinberger et al. (1988, 1989). Changes in the composition and abundance of soil free-living 
nematodes were influenced by the availability of resources (e.g., OM) and adaptation to plant 
eco-physiological traits. A. horridus plants provided litter, food and shelter to the soil animals 
could enrich the biotic components with endoparasites during stopovers of mammals. Moreover, 
we found that soil free-living nematodes, according to the results of 18s rDNA sequences, 
belonged to Brugia, Cooperia, Cylicostephanus, Deladenus, Enterobius and Haemonchus genera 
that were parasites on vertebrates, representing approximately 30% of the taxa, while the 
remaining 70% of the taxa were composed of 37%, 21%, 21% and 21% for BF, FF, PP and OP, 
respectively. These trophic visits by different animals might either be beneficial or cause harm in 
a long-term, due to the fact that animals might be carriers or even trajectories of diseases to other 
animals and even to humans (Jasmer et al., 2003; Holsback et al., 2013). Leung and Koprivnikar 
(2016) found that soil free-living nematodes contributed to the local economical loss as they 
parasitized plants and/or infested a wide range of animals, including humans. 

The result of Venn analysis found that the interaction between both sampling locations and 
plant gender determined the percentage of soil free-living nematodes under each A. horridus plant 
(Fig. 4). Percentages of soil free-living nematodes accounted for 8396 beneath the male, female 
plants and control at Delta, 5096 beneath the female plants at Gobabeb, and 5096 beneath the 
female plants and control at FE and dunes, respectively. According to Venn analysis, percentages 
of soil free-living nematodes at different sites varied: at Delta, 17% species occurred in the union 
of MNFNC; at Gobabeb, 50% species occurred in the unions of MNC, FNC, MO FC; and 5096 
species at FE sand dunes could be explained by the unions of MNF and FNC. This interplay 
elucidated the important effect of plant gender on soil free-living nematodes with the increase in 
aridity in environment. This result was similar to the findings of Beare et al. (1995), Ettema et al. 
(1998) and Yeates and Bongers (1999), who emphasized the importance of spatial scale and plant 
cover in a microhabitat to evaluate the diversity of nematodes. 


5 Conclusions 


The present study showed that soil free-living nematodes were strongly correlated with the 
changes in abiotic components. Moreover, the use of molecular tools had allowed us to increase 
the identification of soil free-living nematodes and endoparasites that were related to mammals in 
such xeric environment. The results from RDA analysis elucidated the importance of soil OM as 
one of the most important variables in determining density and trophic group composition of soil 
nematodes. These results also suggested that some other factors that were not considered in this 
study might contribute to the unexplained variation. 
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